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ABSTRACT 

The power of a relativistic jet depends on the number of leptons and protons carried by the jet 
itself. We have reasons to believe that powerful 7-ray fiat spectrum radio sources emit most 
of their radiation where radiative cooling is severe. This helps to find the minimum number 
of emitting leptons needed to explain the radiation we see. The number of protons is more 
uncertain. If there is one proton per electron, they dominate the jet power, but they could be 
unimportant if the emission is due to electron-positron pairs. In this case the total jet power 
could be much smaller. However, if the 7-ray flux is due to inverse Compton scattering with 
seed photons produced outside the jet, the radiation is anisotropic also in the comoving frame, 
making the jet to recoil. This Compton rocket effect is strong for light, electron-positron jets, 
and negligible for heavy, proton dominated jets. No significant deceleration, required by fast 
superluminal motion, requires a minimum number of protons per lepton, and thus a minimum 
jet power. We apply these ideas to the blazar 3C 454.3, to establish a robust lower limit to 
its total jet power: if the viewing angle 6 V ~ 1/T the jet power is larger than the accretion 
luminosity for any bulk Lorentz factor T. For 6 V = 0°, instead, the minimum jet power 
can be smaller than for T < 25. No more than ~10 pairs per proton are allowed. 



Key words: galaxies: active-galaxies: jets-galaxies: individual: 3C454.3 ■ 
anisms: non-thermal 



radiation mech- 



The high quality data of the the Large Area Telescope (LAT) on- 
board the Fermi satellite, together with the simultaneous observa- 
tions performed by the Swift satellite in the optical-UV and X-ray 
bands and by ground based telescopes allowed a new era in the 
study of blazar jets. Detailed modelling of these sources allows to 
estimate the physical parameters of the jet emitting region, such 
as its magnetic field, particle density, size and bulk Lorentz factor. 
Therefore we can estimate the power that the jet carries in the form 
of particles and fields, and compare it with the accretion luminos- 
ity, at least in Flat Spectrum Radio Quasars (FSRQs) where the disk 
component is visible. In our previous studies of 7-ray loud FSRQs 
(Ghisellini et al. 2010a, 2010b) we found that the jet power can be 
even larger than the accretion luminosity, and it correlates with it 
(also when accounting for the common redshift dependence). 

When estimating the jet power in this way there are two crucial 
uncertainties: i) the total number of leptons, that depends on low en- 
ergy end of the particle distribution (as those are most numerous) 
yet difficult to observe because of synchrotron self-absorption, and 
ii) the number of protons per lepton. For the first concern, evidence 
is accumulating that in FSRQs the radiative cooling is severe, so 
that leptons of almost all energies do cool in one light crossing time, 
and the presence of low energy particles is often required to repro- 
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duce the observed X-ray spectrum, interpreted as inverse Compton 
radiation with photons originating externally to the jet (i.e. External 
Compton, EC hereafter). 

The second concern (how many pairs per proton) has been dis- 
cussed, among others, by Ghisellini et al. (1992); Celotti & Fabian 
(1993); Sikora & Madejski (2000); Celotti & Ghisellini (2008); 
Ghisellini et al. (2010a). If pairs are created in the 7-ray emis- 
sion region we should see a clear break in the spectrum, and the 
absorbed luminosity should be reprocessed at lower energies, es- 
pecially in the X-ray band, where instead the spectral energy dis- 
tribution (SED) of FSRQs has a minimum. If the pairs are created 
very close to the black hole, there is a maximum number of them 
surviving annihilation, corresponding to a local pair scattering op- 
tical depth r± ~ 1 (Ghisellini et al. 1992). When arriving to the 
parsec VLBI scale, the corresponding pair density is less than the 
lepton density required to produce the synchrotron flux we see. On 
the other hand, the 7-ray emitting zone is much smaller and closer 
to the black hole than the VLBI zone, and the number density of the 
surviving pairs might be enough to account for the radiation pro- 
duced in this region. We have found in our earlier works (Celotti 
& Ghisellini 2008; Ghisellini et al. 2010a,b) that the power spent 
by the jet to produce its radiation is often greater than the power 
in Poynting flux and bulk kinetic energy of the emitting leptons, 
requiring an additional form of jet power, i.e. protons. The sim- 
plest hypothesis of one proton per electron leads to jet powers sys- 
tematically larger than the accretion luminosity. So it is crucial to 
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evaluate how many protons per emitting lepton there are in the jet. 
To this end we introduce here a new argument, that can be applied 
when most of the radiation is produced through Compton scatter- 
ing with external radiation (EC). In this case, the emission pattern 
is anisotropic in the comoving frame of the emitting region, that 
must recoil. This is the "hot version" of the Compton drag effect 
(the emitting particles are relativistic in the comoving frame) and is 
called "Compton rocket" (CR hereafter) effect. First studied in the 
'80s, (O'Dell 1981) as a way to accelerate jets using the accretion 
disk photons as seeds, it has been used as a tool to limit the jet bulk 
Lorentz factor T assuming, as seed photons, those re-isotropized 
by the broad line region (BLR) or by a relatively distant torus (e.g. 
Sikora et al. 1996). More recently, it has been used as a decelerat- 
ing agent for very fast jet "spines" moving inside slower jet "layers" 
(Ghisellini et al. 2005), or for large scale jets interacting with the 
cosmic microwave background (Tavecchio et al. 2006). 

We use the CR effect to limit the number of pairs, assuming 
that the jet is moving with a given V, and requiring it does not sig- 
nificantly decelerate by the CR effect, in order to be consistent with 
observations of fast superluminal motion at the VLBI scales. Light 
jets (i.e. pair dominated) can be decelerated more effectively than 
heavier jets (i.e. with an important proton component). Therefore 
requiring no significant deceleration fixes the minimum number of 
protons per lepton, one of the most important number to find out a 
limit on the total jet power. Furthermore, if the (energetically dom- 
inant) 7-ray flux is EC emission, we can in a rather straightfor- 
ward way evaluate the jet powers in its different forms: magnetic, 
leptonic, radiative and protonic, and how these different jet pow- 
ers change by changing V (similarly to what done in Ghisellini & 
Celotti 2001). Doing this, one finds a minimum power, approxi- 
mately where the Poynting flux equals the other dominant form of 
power (i.e. bulk motion of leptons, or protons, or radiative), and a 
corresponding F. 

We apply these arguments to 3C 454.3, one of the best studied 
blazars, used as a test case. It is a FSRQ at z = 0.859 (Jackson & 
Browne 1991), superluminal (Jorstad et al. 2005; Lister et al. 2009) 
with components moving with /3 apP from a few to more than 20, re- 
sulting in estimated bulk Lorentz factors from 10 to 25. It is one of 
the brightest and most variable FSRQs. In April-May 2005 it un- 
derwent outburst, dramatic in optical (Villata et al. 2006) and vis- 
ible also at X-ray energies (Pian et al. 2006; Giommi et al. 2006). 
The AGILE satellite detected 3C 454.3 as one of the brightest and 
variable sources in the 7-ray band (Vercellone et al. 2007, 2010). 
After the launch of the Fermi satellite 3C 454.3 was seen to flare 
several times (Tosti et al. 2008; Abdo et al. 2009), with a climax in 
December 2009 (Bonnoli et al. 2010). 



2 THE COMPTON ROCKET EFFECT AND THE PAIR 
CONTENT OF THE JET 

If the main emission process of blazars is synchrotron and self- 
Compton radiation, then the emitted luminosity is isotropic in the 
comoving frame. This means that the jet loses mass, but not ve- 
locity. The lost mass is at the expenses of the "relativistic" (i.e. 
7m.eC 2 ) mass of the emitting leptons. Instead, if powerful blazars 
produce most of their emission by scattering external radiation, the 
produced radiation is anisotropic even in the comoving frame, and 
the jet must decelerate (i.e. it recoils in the comoving frame). The 
amount of this deceleration depends on the produced EC luminos- 
ity and by the inertia of the jet, i.e. if the jet is "heavy" or "light". 
Therefore this "Compton rocket" (CR) effect (i.e. the "hot" ver- 



sion of Compton drag, because the leptons are relativistic) can give 
some limit on the minimum number of protons present in the jet, 
requiring that is does not decelerate significantly. We here present a 
simple derivation of the relevant formulae, that agree with the more 
detailed and complex derivation of Sikora et al. (1996). 

Let us remain in the observer frame. There we measure an ex- 
ternal and isotropic radiation energy density £/ cxt ■ The total Lorentz 
factor (7) of the electrons is the superposition of the bulk (r) and 
random (7) Lorentz factors (/Sbuik and j3 are the corresponding ve- 
locities). Assume that the bulk motions occurs along the x-axis, 
and that the random velocity forms an angle 8' with respect to that 
axis, in the comoving frame. We have (e.g. Rybicki & Lightman 
1979): 



P' cos 9' + /3 bu]k _ 
l + /3buik/3' cosfl'' 

The total 7 2 is 

7 2 = 



y r(l+/3 bulk /3'cos6") 
(l + /3buik/3' cos0')Vr 2 



(1) 



(2) 



If the particle distribution is isotropic in the comoving frame, the 
average over angles gives 
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which gives the factor (4/3) for ultra-relativistic speeds. 

Now assume that a portion of the jet carries a total number N p 
of protons and iV c leptons (including pairs). The "cooling time" of 
the jet (i.e. the time for halving T) is 



^cool 



E N p m p c 2 + N c {-y)m a c 2 

E ~ (4/3)<TTCATet/ext(7 2 ) 

9 (N p /N c )m p c 2 + (7) m c c 
16 ' 



CTTCt/ext<7 2 )r 



(4) 



In the cooling time t C ooi, the jet travels a distance R coo i = ficbmol- 
The corresponding interval of time as measured by the observer 
is Doppler contracted by the factor (1 — /3cos# v ) = 1/(F5) 
(6 V is the viewing angle and 8 the beaming factor). The time 
tcooi(l —fi cos V ) has to be compared with the timescale for which 
the particles are indeed described by an energy distribution with the 
value of (7 2 ) used, provided that, during this time, the radiation 
density remains [/ cxt • This timescale is approximately the variabil- 
ity timescale. The CR effect is unimportant if 
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(l + z)(l- P cos 6 V ) (1 + z) 
9(1 + z) (N p /N c )m p c 2 + (7)m c c 2 



16 
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(5) 



This limit becomes very severe if the jet is dominated by hot pairs 
and if Uext = (7blr, the radiation energy density is dominated by 
radiation from the broad line region. In this case jets with T ^ 10 
are bound to decelerate. They do not decelerate if they contain a 
proton component that increases their inertia. We can rewrite Eq. 
[5] to find the minimum ratio N p /N e compatible with halving V in 
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OTct/ cxt <7 2 } (7)? 



9 1 + z 



(6) 



Since N c = N± + N p (pairs plus electrons associated with pro- 
tons), we have N±/N p = (N a /N p ) - 1. 
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3 THE POWER OF THE JET 

The most robust estimate on the jet power is the power P r spent 
to produce the radiation as measured with a "47r" detector sur- 
rounding the source in the observer frame. If the luminosity 
L' is isotropic in the comoving frame we would derive P r — 
L7(4tt) JT S 4 dfl = (4/3)r 2 L'. But if the EC process is impor- 
tant the emission is not isotropic in the rest frame, and the observed 
flux, instead of being boosted by <5 4 , follows a pattern given by 
5 i (6/T) 2 (see Dermer 1995 and Georganopoulos et al. 2001). Set- 
ting (L') the angle averaged luminosity in the comoving frame, we 
then have 

(L>) f S 6 (6) 16 2 16r 4 L obs 

The power in bulk motion of leptons, protons, and magnetic 
fields are calculated as: 



nRbiobT 2 (3c m e c 2 
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where Pbiob is the size of the emitting source. We are assuming 
that protons are cold and that n c is the total number density of 
leptons, including pairs (if present), so that n c = n± + n p . We 
also assume that all leptons are relativistic and are described by the 
energy distribution N(j). Neglecting cold leptons minimizes the 
power requirement. The total jet power is Pj e t = Pr+Pc+Pp+Ps- 
As long as the scattering is in the Thomson regime the ob- 
served luminosity in the EC component of the SED is: 

l ec ~ g o- T cn c (7 )U cxt j^S (9) 

We can then find the number density n c of the emitting leptons: 



9L E c 



16^Lb^Tc{7 2 }t/ext«5 6 



(10) 



The averages (7) and (■y 2 ) are calculated assuming the emitting 
particle distribution is a broken power law, extending from 7 coo i to 
7 pca k with slope TV (7) oc 7~ 2 , as appropriate for radiative cooling, 
and breaking above 7 poa k, where we assume TV (7) oc 7 _p up to 
7max- The slope p is related to the observed energy spectral index 
q above the synchrotron and the EC peaks as p — 2a + 1. The 
values of 7 CO oi and 7 pea k will depend on T and S (see Eq.ll5land 
Eq.|16l>. In general, both {7) and (7 2 ) decrease by increasing T and 
5: this is because the external photon field is seen more boosted in 
the comoving frame, inducing a stronger Compton cooling (and so 
7cooi decreases); at the same time a smaller 7 pea k is required to 
produce the high energy peak. 

The magnetic energy density Ub can be derived in terms of 
the "Compton dominance", namely the 7-ray to synchrotron lumi- 
nosity ratio L 7 /L 8yn : 
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4 THE MINIMUM JET POWER 

We now show how the different forms of jet power changes by 
changing the bulk Lorentz factor. These estimates depends on the 



following parameters: L 7 ~ Lec; L 1 /L sy - n \ v c (the vF v peak 
frequency of the 7-ray spectrum); V (the viewing angle); t val — 
(l+z)Rbiob/{c5); a~f (the energy spectral index of the 7-ray spec- 
trum above the peak) and z. Of these 7 parameters, all but one 
(the viewing angle) are observables. Note that the external radia- 
tion energy density is not a free parameter if the typical radius of 
the BLR (or the reprocessing torus) depends on the disk luminosity 
as _R cxt oc L X J 2 . In this case [Text oc Ld/R 2 xt I s constant. 

The viewing angle V is in general unknown. We can assume 
8 V = 1 /r, in order to always have V = 5. Alternatively, we can 
assume 9 V = 0, i.e. S — 2T. Maximizing the Doppler boosting, 
this choice will minimize the derived powers but the probability to 
observe any source with 6 V — is vanishingly small. There might 
be an exception: consider the case for which the velocity vectors of 
the emitting flow are not parallel, but somewhat radial within the 
jet aperture angle #j e t. Thus if 9 V is finite, but smaller than 6j e t, 
there is a portion of the jet exactly pointing at us. On the other 
hand, the corresponding emitting volume of this portion is small, 
and the flux we see is mostly contributed for by those parts of the 
jet moving with 9 V ~ 1/T, because they have a larger volume. 

We are led to conclude that the case 6 V = 1/T is favored. It 
is also the angle for which the superluminal motion is maximized, 
but this implies that the jet does not change direction between the 
7-ray emitting region and the VLBI scale, that is not guaranteed. 

We now rewrite the different forms of jet power in order to 
make more transparent their dependences on the bulk Lorentz fac- 
tor and the parameters listed above. 
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i p from Eq. [6] 



(12) 



The last approximate equality assumes that the second term in 
round brackets in Eq.[6]is negligible with respect to the first. Note 
that, as expected, the minimum P p limited by the CR effect is of 
the order of the power spent in radiation (given in Eq.[7]>- B e aware 
that (7) and (7 2 ) do depend on F (see Eq.[l5]and Eq.ll6lbelow). 

4.1 Application to the blazar 3C 454.3 

In order to find a lower limit to the jet power of 3C 454.3 we make 
the following assumptions: 

• In Bonnoli et al. (2010) we have shown that during the big 
flare of Nov-Dec 2009 the 7-ray, X-ray and optical fluxes of 3C 
454.3 were correlated with one another, with the 7-ray flux varying 
more than linearly with the flux in the other two bands. We take 
this as a very robust indication that most of the non-thermal flux 
received from 3C 454.3 above the far IR band is produced in the 
same region of the jet. 

• The size of the emitting region i?biob is assumed to be asso- 
ciated with the minimum variability timescale t var of the source. 
In the 7-rays, Tavecchio et al. (2010); Foschini et al. (2010) and 
Ackermann et al. (2010) found significant variations in 3-6 hours. 
Therefore 

R ^=«™ih~ 7xi ° 1B (^)Q cm (13) 

• We assume an accretion disk luminosity Ld ~ 6.7 x 10 46 erg 
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Figure 1. Different form of jet power as a function of the bulk Lorentz factor 
r, as labelled. For this particular example, we have assumed L 7 = 9 X 10 49 
erg s , v c = 10 22 Hz, t var = 6 hours and L 7 /L syn = 20, appropriate 
for the blazar 3C 454.3, at z = 0.859. We have further assumed that 6 = I\ 
implying 9 V = 1/T. The horizontal grey line indicates the accretion disk 
luminosity. The long dashed lines correspond to Pj c t assuming one proton 
per emitting electrons (i.e. no pairs) or instead assuming the minimum num- 
ber of protons per electron consistent with no strong jet deceleration for the 
Compton rocket effect. We can see that Pj c t > Ld for all T. The bottom 
panel shows the maximum pair to proton ratio allowed by the CR effect. 



s _1 , based on direct detection of the Lyman-a line (Bonnoli et al. 
2010) and on the flattening of the optical-UV SED when the source 
is in low state. 

• We assume that the BLR reprocesses 10% of Ld and that the 
BLR size is given by Pblr = 10 17 Lj % 5 cm. This choice (in rough 
agreement with Bentz et al. 2006 and Kaspi et al. 2007), implies 
that the radiation energy density within the broad line region is con- 
stant: 

O.ILd 1 



C^BLR = 



4ttP 2 



12tt 



erg cm 



(14) 



The short t var suggests that dissipation takes place within the BLR, 
so we assume (7 cxt = (7blr- 

• After one light crossing time Pbiob/c = t V ar<5/(l + z) the 
cooling energy 7 coo i is 



7cool 



3(1 + z)m e c 2 



4.0-TCivar 5 [U^ hR + U B + f/s'yn] T 2 <5 



(15) 



with the EC mechanism (with BLR photons as seeds) being the 
dominant cooling agent. 

• The soft slope of the 7-ray spectrum and the hard slope of the 
X-ray spectrum constrain the peak of the high energy component 
hu c of the SED to lie close to 100 MeV. For the EC process, the 
peak is made by electrons at 7 pea k scattering the Ly-cv seed pho- 
tons with frequency vl V c. Then 7 pGa k is given by 



Figure 2. Same as Fig.[T] but assuming 8 V = 0. Note the different range 
of the y axis. In the bottom panel, the function n±/n p is similar, but not 
identical to the one shown in Fig.fTJ 



r<5 
T+z 



Tpeak 



M 1 + z ) 



1/2 



(16) 



For all reasonable parameters appropriate for 3C 454.3, 7 pca k > 
7cooi, implying that most of the energy injected in the form of rel- 
ativistic leptons is radiated away in one light crossing time. 

The top panel of Fig. [TJ shows the different jet powers as a 
function of V assuming that 9 V — 1 /F. The two dashed lines are the 
total Pjct derived assuming one proton per electron (i.e. no pairs) 
or instead the minimum number of protons given by Eq. [6] In both 
cases Pj ct > Ld for all values of F. The minimum Pj c t is set by 
the equipartition between P p and Pb. This occurs at F ~ 55 for the 
"no pairs" case, and F ~ 40 for the case of the minimum number 
of protons ("no decel." case). 

The bottom panel of Fig. [TJshows the maximum ratio n± /n p 
required not to decelerate significantly for the CR effect, as a func- 
tion of r. To understand the behavior of this curve consider Eq.[6] 
If we neglect the second term, we have n±/n p oc (r 2 <5{7 2 )) -1 . 
For illustration, consider a particle distribution extending only be- 
tween 7cooi and 7 pca k with slope ^(7) oc 7~ 2 . In this case 
(7 2 > ~ 7cooi7 P cak oc [r 2 5(r<5) 1/2 ]~ 1 as can be seen through Eq. 
I15land Eq.[T6] Therefore n±/n p oc (F3) 1 / 2 oc F. This behavior 
ends when 7 coo i becomes unity, i.e. for large values of F. In this 
case n±/n p oc F~ 1 (FS)~ 1 ^ 2 ~ F~ 2 . The maximum in n±/n p 
therefore occurs when 7 CO oi becomes unity. 

The top panel of Fig. [2] shows the powers assuming 9 V — 0° . 
In this case the minimum Pjct occurs for r ~ 18 (no pairs) or 
r ~ 16 (with N p /N c given by Eq. [6}. In the case of no pairs 
Pjct ~ Ld, and is a factor 4 smaller than Ld for the minimum 
number of protons allowed by Eq. [6] The ratio n± /n p behaves 
approximately as in Fig.Q] 
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5 CONCLUSIONS 

It is very likely that the 7-ray emission region in powerful FS- 
RQs is within their broad line region, with broad line photons be- 
ing the seeds for the inverse Compton scattering process. We base 
this assumption on the observed fast variability, difficult to explain 
in models where dissipation takes place at much larger distances in 
the jet, as in the models by Marscher et al. (2008) and Sikora et al. 
(2008; 2009). This implies that the radiation is anisotropic in the 
comoving frame, making the jet to recoil. The observer would then 
see a deceleration of the jet, important for pure electron-positron 
light jets and becoming less significant if the jet is heavier due to 
the presence of protons. Therefore the requirement of no or only 
modest deceleration translates in a requirement on the amount of 
protons in the jet. This then gives a lower limit on the total jet 
power. 

Within the framework of synchrotron and external Compton 
models, the only parameters that remain somewhat free (namely 
not accurately given by observational data) for calculating the min- 
imum Pj t are the viewing angle 8 V and the bulk Lorentz factor 
r. We can however see how the minimum Pj at values change as a 
function of T, assuming a given viewing angle. Doing so, we find 
at which T the jet power is minimized. We can then compare this 
minimum P™ t m ("minimum of the minimum values") with the ac- 
cretion disk luminosity. 

Applying these arguments to 3C 454.3, one of the best studied 
7-ray blazars, we found that if 8 V — 1/T, then P™ t m > Ld, while 
it becomes a factor 4 smaller than Ld if 6 V = 0° . The key ques- 
tion if the jet power can be larger than the accretion disk luminosity 
remains therefore open, but with a narrower range of possibilities 
than before. We can exclude pure electron-positron jets, but we can 
allow for ~ 10 pairs per proton. This value is in agreement with 
what found by Sikora & Madejski (2000) using a different argu- 
ment. Be aware that all these estimates are based on the assumption 
that all leptons present in the source participate to the emission. If 
cold leptons were present, they would increase our estimate of the 
jet power. 
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